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ABSTRACT 


A "typical eddy" box model for the second-order closure of 
turbulent, multi-species, reacting flows has been developed. 

The model structure is quite general and is valid for an arbitrary 
number of species. For the case of a reaction involving three 
species, the nine model parameters are determined from equations 
for the nine independent first- and second-order correlations: 

a , 6 , T , a 1 g ' , a'y' , B'y' , a"' T 1 , I 7 ! 71 " , and T'T* . _ The 

model enables calculation of any higher-order correlation involving 
mass fractions, temperatures, and reaction rates in terms of first- 
and second-order correlations. Model predictions for the reaction 
rate are in very good agreement with exact solutions of the reaction 
rate equations for a number of assumed flow distributions. 


INTRODUCTION 


Recent advances in the analysis of turbulent flows by closure 
of the transport equations for various Reynolds stresses at the 
second-order level have resulted in a deeper understanding of the 
physics of turbulent flows that could not be obtained by classical 
first-order closure methods such as the eddy viscosity method. 
A.R.A.P. has been in the forefront of the application of second- 
order modeling techniques to the calculation of turbulent shear 
flows. We have successfully applied these techniques to a variety 
of flow problems including incompressible and compressible free 
shear layers and boundary layers (ref. 1), simulation of the flow 
in a HF chemical laser (ref. 2), behavior of the planetary boundary 
layer (ref. 3), and the decay of axisymmetric vortices (ref. 4). 

In the course of these studies, we have developed and tested models 
for a large number of third-order correlations that appear in the 
transport equations for the second-order correlations. These models 
will have direct application in our current development of a 
computer program to analyze turbulent reacting flows. 

The analysis of turbulent, multi-species, reacting flows is an 
example of the kind of problem that requires a more sophisticated 
approach than the eddy viscosity method. During the past two years, 
A.R.A.P. has been developing an invariant second-order coupled 



diffusion and chemistry model with support from NASA, EPA, and DOT. 
The model has been used to study a number of atmospheric (incom- 
pressible) chemistry problems such as the dispersal of pollutants 
from industrial stacks and highways, the behavior of the far wake 
of the SST exhaust, and other similar problems (ref. 5)- These 
studies demonstrated that there are indeed real world atmospheric 
pollution problems in which neglect of the fluctuations of 
concentrations of reacting species introduces significant errors .■ 
Neglect of the effects of inhomogeneous mixing could result in 
order of magnitude errors in the chemical reaction rates. An 
approximate closure scheme for the chemical submodel which conforms 
to the principles of invariant modeling and which accounts for the 
effects of inhomogeneous mixing over a wide range of conditions was 
developed. 

The above-mentioned studies of modeling of reacting flows were 
limited to problems involving low heat release in which the fluctua- 
tions in the reaction rate could be neglected. However, the success 
of these studies raised the distinct possibility that by combining 
the compressible flow models with further advances in chemical 
modeling, it will be possible to construct a second-order closure 
model for the computation of compressible flow problems involving 
highly exothermic chemical reactions, such as the analysis of a 
turbulent diffusion flame. A.R.A.P. is currently being funded by 
NASA under Contract NAS1-12412 to develop the models for third-order 
correlations involving fluctuations of species concentrations , 
temperatures, and reaction rates. This report details the develop- 
ment of these models over the past year. 


CHEMICAL KINETIC EQUATIONS 


Consider an irreversible single-step reaction between chemi- 
cal species a and 6 to form y . 

a + 6 -* y ( 1 ) 

We will primarily deal with such a three-species system, although 
the model can easily be extended to reactions involving more 
species; for example, 

a + B -*• y + 6 

The reaction rate equations may be written 

= -k[a ] [ B ] 

( 2 ) 

= -k[a ] [ B] 


2 



If, instead of the molar concentrations [a] and [3] , 

Eqs . (2) are written in terms of the mass fractions a and 3, 
the reaction equations become 

5tT ' k i“ 8 k i = k 

ft = - k 2° 6 k 2 * k iT <4) 

a 

where p is the density, W and W D are the molecular weights 

Ot p 

of species a and 3 , respectively, and k-^ and k^ are the 
reaction rate constants. d 

Equations (3) and (4) specify the local instantaneous rate 
of change of the concentration of the reactants. In a turbulent 
reacting flow, these equations may be written in terms of the mean 
and fluctuating values by writing a given quantity <j> as the sum 
of a mean and fluctuating part 

<J> = $ + <j>' (5) 


If Eqs. (3) and (4) are expanded by the scheme given in (5), the 
result is 


3a 3a* 


3t 3 1 


-(^ + kp(a + a')(3 + 3') 


= ~(k 1 + k^)(a3 + a3' + 3a' + a'3’) 


( 6 ) 


and 


36 , 33' 

3t + at " 


-(k 2 + kp(a3 + a3' + 3a' + a'3') 


(7) 


The chemical kinetic equations for the average rates of change of 
the concentrations of a and 3 then become 


and 


3a 
3 1 


k^(a3 + a'3') - a k^3' 


3 k^a • - k^a»3’ 


( 8 ) 


33 
3 1 


k^CotS + a*3') - a k^3* 


3 k^a 1 - k^a * 3 ' 


(9) 
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In a second-order closure scheme, one generally carries 
equations for all the independent second-order correlations of 
the fluctuations. Chemical kine tic equations can be easily 
derived for terms such as a 'B' as follows. By subtracting 
Eqs . (8) and (9) from (6) and (7), respectively, one obtains 

|^- = - k 1 a6’ - k 1 Ba' - k-^a’B’ - a' 3' ) - aBkj^ 

- a ( k ' 3 ' - k'B' ) - B(k»a» - k'a' ) - (k'a'3' - k'a'3') 


- k^aB' - k^Ba' - k2(a'3' - a'B' 


- aBk^ 


- a(k^3 ' - ki,B ' ) - B(k£a' - k£a') - (kJci'B' - k£a'3') 

(11 

Multiplying (10) by 3' , (11) by a' , and averaging the sum 

of the resulting expressions 


3 — nr 0 , sot* . , 3 3 ' 

at a 6 " 6 at a at 


k 1 (a3’ 3' + 3a’ 3' + a’3' 3’) 


aBk'B' - a k^B’B' - 3 kja'3' - k’a'B'B' 


k^aa'B' + 3a' a* + a'a'B') 


aBk^a' - ak^a'B' - Bk^a'a' - k^a'a'B' 


Expressions for other independent second-order correlations 
can be obtained in a similar fashion. Equations (8), (9), and 
(12) contain higher-order correlations, such as the terms 

a'B'B* , k^a’B* , k^a'B' , k^a'a' , k^a'B'B' , etc'., which now 

have to be modeled in terms of the mean quan tit ies a , 3 , T , 
etc . , and the second-order correlations a'B' , a'a* , B'Y 1 > 
a'T' , etc. Once all the higher-order terms are modeled, we 
have a closed system of equations for the means and the inde- 
pendent second-order correlations which can be solved. We now 
have to develop suitable models for the various higher-order 
correlations . 



"TYPICAL EDDY" BOX MODEL 


In the second-order closure of the chemical submodel for 
a reacting mixture of three species, let us identify what 
information is available to us to construct a general model for 
the higher-order correlations and what the constraints are on 
the model. 

Some effort was devoted to the development of analytical 
model expressions for the third-order correlations following 
procedures similar to those used in reference 5 for low heat 
release reactions. These were not entirely successful b ecause 
of the lack of proper bounds on third-order terms like k’a'a’ . 
It was decided to construct a "typical eddy" for modeling 
concentration, temperature, and reaction rate fluctuations. A 
similar concept is used for modeling other correlations. For 
example, dissipation of turbulent kinetic energy is modeled as 
though all the dissipation was due to eddies of a single scale 
length. The tendency towards isotropy is modeled in terms of 
another scale length, and so on. 

First, let us consider how many variables are at our disposal 
to construct the model of a "typical eddy" for the chemistry sub- 
model. We have the equations for the mean terms 

a , 3 , and T (13) 

Note that y is not an independent variable since 

a + 3 + y = 1 

Further, we have equations for 


a’ 3' , a ' y ’ , 3'Y* , a’T* , 3'T’ , and T' 2 (14) 


These are the only species and/or temperature second-order 
correlations that are independent since a* + 3' + Y' = 0 
and, therefore. 


, f 2 

= - (a ' 3 ' 

+ 

a ' y ' ) 

1 1 2 

= - (a ' 3 ' 

+ 

b’y' ) 

, t 2 

= - (a ' y ' 

+ 

3 ' Y ' ) 

T' 

II 

1 

P 

+ 

3'T' ) 

we 

have a'T 

t 

or J 7 


since the reaction 
rate terms kj_ and k 2 are known functions of T (for example, 
an Arrhenius relationship), we should be able to write expres- 


sions for k^a' , 


k|3' , etc. One of our problems in selecting 


5 


3 


etc. in 


a model is to obtain expressions for k^a* 


kja'B' , 


terms of the known quantities (13) an d (1 4) whi ch are completely 
consistent with our expressions for a'T’ and a’B'T* . 


As listed above, we have nine independent parameters (Eqs. 

( 13 ) and (14)) to construct a typical eddy for a mixture of. 
three species. If we remain at a particular place in a fluid 
and watch the turbulent flow go by and ask, "What does a typical 
eddy look like?" or "What is the typical repetitive structure of 
the flow?", the answer will be some variation of the following. 
Consider two unmixed streams of a and 3, which mix and react 
to form y . In the early stages of the mixing, one will observe 
large amounts of pure a and 3 . Some reaction would have 
occurred and some pure product y may pass by. There would be 
diffusion taking place and we expect to observe fluid elements 
consisting of differing proportions of a + 3 and a + 3 + Y 
(wherein the reaction is proceeding), a + y and 3 + Y • At a. 
later time, one would expect less frequent fluid masses consisting 
of pure a and 3 but the basic structure of the eddy will be the 
same. With this physical and intuitive idea of the flow, we now 
try to construct a model. The degrees of freedom of the model are 
limited by the number of independent parameters available. 


Let us also set down other desirable features of the model 
and certain constraints the model must satisfy. 

1. The three-species model should be capable of 
collapsing to the case of one or two species 
and be logically extendible to larger number 
of species . 

2. The model should not have to label the three 
species; that is, which are the reactants and 
which are the products? Further, the model 
should be capable of handling a nonreacting 
mixture of three constituents. 

3. The model must satisfy the_end of the reaction 
constraints; that is, if a and 3 are the 
reactants , 


a ' 3 ' = - a • 3 


a • 2 3 ' = a 2 3 - a ,2 B 

Consider now the model for a mixture of three species. In 
the absence of any temperature fluctuations, there are only five 
independent parameters with which to construct the typical eddy. 
These are 


a , 3 , a'B' , a ' y ' , and 3 1 Y ’ (15) 
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A desirable "typical eddy" structure, shown in Figure la, 
requires as many as 11 parameters, e g and k^ -* kg . 

Obviously, we cannot have this general a model. A large number 
of simpler models were investigated. The proportions of the 
various materials in cells 4 through 7 can be fixed first. If 
cell 7, in which a, 6 , and y exist together, is part of the 
model, cells 4, 5, and 6 could be considered to have either 
equal amounts of the two constituents or have them in proportion 
to the mean concentrations. In order to be able to consider the 
case of nonreacting mixing of three species for which the end 
result is a uniformly mixed flow, the concentrations in cell 7 
must be proportional to the mean values. If cell 7 is taken 
out of the model (this cell is only absolutely necessary for 
accurate computation when the reaction rates are slow compared 
to diffusion processes), the concentration in cells 4, 5, and 6 
must be in proportion to the mean values, for the model should 
be capable of handling a confined flame in which the reactants 
are not available in stoichiometric proportions. Then A the end 
result is a uniformly mixed flow of a and y (or $ and y ) 
and the amounts of a and y in cell 5 should be proportional 
to the means. By symmetry, we must do the same for cells 4 and 6. 

With the above-noted assumptions, we have six unknowns, 

-* Eg , in the model. Various possible models were investigated. 

Some of the models which were rejected are sketched in Figures lb 
and 1c. These failed for different reasons — lack of symmetry, 
not satisfying the reaction end limit, etc. 

After extensive studies, we were left with two models which 
met all the conditions and tests we applied to them and appeared 
to be quite satisfactory. The predictions of these models were 
compared to the true chemistry results for a number of assumed 
flow distributions, and the results were very good. The models 
and the tests are discussed later in this report. 

Both of these models were a significant improvement over 
"classical" chemistry calculations, but in some flows Model I 
showed better results while Model II did better in other flows. 

A seven cell, six parameter model combining the better features 
of these two models appeared to be desirable, but it would 
require an additional independent parameter, and we could not 
find a consistent way of introducing another parameter. In the 
last month we have developed a new seven-cell model (Model III) 
which still requires only five parameters. This model is 
discussed in a later section "Description of Model III." It must 
be tested more completely in the same way that Models I and II 
were tested, but the final model is so closely related to these 
two models that we expect it to show very similar results. 
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DISCUSSION OF MODELS I AND II 


We now consider the joint species temperature models. 

With no temperature fluctuations, there are five- parameters 

(a, 3 , a’ 3' , a'y* , 3'y' ) for the three-component species 

model. For variable temperatures, we have four additional 

parameters (T , T ' T ' , a "I ' , 3'T' ). For a mixture of two 

species, we have a total of five parameters: a , a'B* , T 

T'T' # a ' T ’ . For the case when only one species is present, 

the model must admit temperature fluctuations whose distribu- 
tion functions are almost Gaussian, and to allow this result we 
assume that in the temperature distribution functions each cell 
spends one-half its time at a temperature AT above its base 
temperature. In the case of one species, the AT is governed 
by wall boundary conditions and/or viscous heating within the 
shear layer. The one- and two-species models are shown in 
Figure 2. 

The six cell, three-species models I and II are shown in 
Figures 3 and 4. The temperature model has four unknowns: T a , 

T^ , T , and AT . The temperatures of the other three cells 

are the appropriately weighted means of the three base temper- 
atures for the three species. 

The concentration model and the temperature model allows 
one to compute any higher-order correlation of scalar fluctua- 
tions involving concentrations or temperatures. The second and 
higher-order correlations involving the fluctuations of the 
reaction rate constant k' can also be calculated, because it 
is possible to calculate the instantaneous reaction rate k and 
plot the reaction rate distributions for each half ce.l l a s sh own 
in the figure. Then it is possible to calculate k'a' , k'3' or 
any other higher-order correlation of k' with fluctuations in 
concentration or temperature. 

The three-species models described above have been compared 
against exact solutions of the reaction rate equations. The 
procedure is as follows. A number of base flows (assumed concen- 
tration and temperature distributions) are set up. A computer 
program has been written to calculate the nine model parameters 
from the values of the nine first- and second-order correlations 
for a given base flow at specified modeling times. The program 
then calculates various moments for the base flow as the reactions 
proceed in time. A comparison of the model predictions and the 
base flow (or true) chemistry is displayed using a plotting 
routine . 

The following reaction rate expression has been used for 
the tests reported here: 
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^p 5 "^*5 = a 5 Cpa"*" a+ ^5^p y~^ Y 

Cp 6 T 6 = ^ 6 C p ^T^+ x 6 C Py T y 



note K # k ( T ) 


0 £i fl 

2 2 

Figure 4. Model II for a typical eddy in three-species flow 
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ic n 1 4 ( 1.600 \ 

= 1-5 x 10 exp ^ WF~ ) 

This corresponds to the reaction 

a + g -*■ y 


H,. + F ■* HF* + H 

d 

in a HF chemical laser, A.R.A.P. is studying the turbulent 
mixing and reaction in these devices with support from ARPA . 

Figure 5 shows first the influence of various terms in 
the expression for 3a/ 3 t (eq. (8)). An exact solution of 
the chemical kinetic equations was obtained for base flow 1, 
and all tne moments of interest were generated. The graph 
clearly illustrates the importance of taking the second-order 

k'a' and k*8' , and the third-order k ' a ' B ' correlation 

into account for predicting the reaction rate. 

Figures 6 and 7 present the results of modeling base flow 
1 for the variable temperature and the constant temperature 
cases. The parameters of the base flow are shown on the figure. 
The "classical" chemistry results are obtained as 3a./ 3 1 = 
kCTj * a ■ 8 ; that is, by ignoring the turbulent fluctuations. 
Models 1 and II are established at time zero from the parameters 
of the base flow distribution function and the models are then 
allowed to proceed in time. It can be seen that the inclusion 
of temperature fluctuations leads to larger differences between 
the classical and true results. The model predictions in both 
cases show a very significant improvement over the classical 
results, but Model I is definitely superior for this flow. 

Results for three other base flows are presented in Figures 
8 through 1 3 . Fop base flow 2, the selected temperatures are 
quite high, and the reaction rate values are virtually saturated. 
Consequently, Figures 8 (variable temperatures) and 9 (constant 
temperatures ) do not exhibit large differences. The model flows 
are markedly better than classical chemistry calculations in 
both cases , but Model I does better In one case 'while Model II 
follows better the true chemistry result in the other case . For 
base flow 3, the results shown in Figures 10 and 11 indicate that 
Model I is superior to Model II. • Finally, tor base flow 4, we 
see that for the variable temperature case (fig. 12), Model II 
would be the preferred model, but the constant temperature 
studies .in Figure 13 suggest that Model I is the better model. 

In all. these cases studied, the model predictions are 
significant ly better than classical chemistry calculations and 
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Figure 5. Component terms in the expression for rate of change of 



Figure 6. Base flow 1 to check correlation modeling; variable 
temperature cells 
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Figure 7. 


Base flow 1 to check correlation modeling; constant 
temperature cells 
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Time, t 

Base flow 2 to check correlation modeling; variable 
temperature cells 






Figure 9 - 
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Base flow 2 to check correlation modeling; constant 
temperature cells 
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Figure 10. Base flow 3 to check correlation modeling; variable 
temperature cells 
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True chemistry, 
Model I 







Figure 11. Base flow 3 to check correlation modeling; constant 
temperature cells 
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Figure 13 
22 


Base flow 4 to check correlation modeling; constant 
temperature cells 
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one or the other of the two models shows good agreement with 
the exact solution of the chemical kinetic equations. We have 
now developed a seven-cell model which incorporates the better 
features of the two models discussed above. This model. Model 
III, is described in the next section. 


DESCRIPTION OP MODEL III 


For this model, the typical eddy of extent unity in time 
is considered to have the following structure. For a time eq 
the eddy contains only the species a . For a time £5 the 
eddy contains only the species 6 , and for time £3 the eddy 
contains only the species y . For time equal to £4 , £5 , and 
eg , the eddy is assumed to contain 1/2 a and 1/2 8 , 1/2 a 
and 1/2 Y , and 1/2 6 and 1/2 y by mass, respectively, in a 
state of molecular mixedness. In addition, we assume that for a 
time 



the eddy contains a , 8 and y uniformly mixed and in propor- 
tion to the mean mass fractions a , 8 , and y . Since the 
total extent of the eddy is unity. 


£ 1 + "2 + e 3 + % + e 5 + e 6 = (1 " e 7 




g ' Y ' | 

By I 


Thus, the model has only five unknowns necessary to define the 
eddy which are determined by the five available equations for 

a ,8 ,a'8' , a' y' , and B'y' . The temperature and the 

reaction rate distributions are set up in the same way as for 
Models I and II. Figure 14 illustrates the complete seven-cell 
model . 


This model is currently being tested by comparison of model 
predictions to exact solutions of the chemical rate equations. 
Basically it is a combination of the two previous models, and we 
expect the results of these tests to be quite satisfactory. The 
model can be extended to the case of four species in a straight- 
forward manner and collapses correctly to the cases of flows 
involving one or two species. 
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Cp 6 T6 s ^eC p ^T^+ y 6 CpyTy 
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Figure 14. Seven-cell, three-species "typical eddy" model 

MODEL III 


00 . 1100 . 



CONCLUSIONS 


A "typical eddy" box model for the second-order closure 
of turbulent reacting flows has been developed. The model 
structure is quite general and is valid for an arbitrary number 
of species. Model predictions are significantly better than 
classical chemistry (involving only the means) computations and 
are in good agreement with exact solutions of the reaction rate 
equations for a number of assumed flow distributions. 

The model enables calculation of any correlation of scalar 
fluctuations involving concentrations, temperatures, and 
reaction rates. It provides the capability of making a second- 
order closure computation of turbulent flames and other flows 
involving highly exothermic chemical reactions. 
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